Mon. Not. R. Astron. Soc. 000, hHia (1999) Printed February 1, 2008 (MN WI^ style file vl.4) 



The effect of a cusped matter distribution on the 
formation of brightest cluster members 



E. Athanassoula^, A. Garijo^ and C. Garcia Gomez^ 

^ Observatoire de Marseille, 2 place Le Verrier, F-13248 Marseille Cedex 4, France 
^ Dep. Enginyeria Informdtica i Matematiques, Escola Tecnica Superior d'Enginyeria, 
Universitat Rovira i Virgili, 43006 Tarragona, Spain 



o 
o 
o 

(N 

Oh. 

l> 



> 

o 
o 
o 

o 



1 



Accepted . Received 



ABSTRACT 

We present N-body simulations of galaxy groups embedded in a common halo of 
matter. We study the influence of the different initial conditions upon the evolution 
of the group and show that denser configurations evolve faster, as expected. We then 
concentrate on the influence of the initial radial density profile of the common halo 
and of the galaxy distribution. We select two kinds of density distributions, a singular 
profile (modeled by a Hernquist distribution) and a profile with a flat core (modeled by 
a Plummer sphere). In all cases we witness the formation of a central massive object 
due to mergings of individual galaxies and to accretion of stripped material, but both 
its formation history and its properties depend heavily on the initial distribution. In 
Hernquist models the formation is due to a "burst" of mergings in the inner parts, 
due to the large initial concentration of galaxies in the center. The merging rate is 
much slower in the initial phases of the evolution of a Plummer distribution, where 
the contribution of accretion to the formation of the central object is much more 
important. The central objects formed within Plummer distributions have projected 
density proflles which are not in agreement with the radial profiles of observed brightest 
cluster members, unless the percentage of mass in the common halo is small. On the 
contrary the central object formed in initially cusped models has projected radial 
profiles in very good agreement with those of brightest cluster members, sometimes 
also showing luminosity excess over the r^' ** law in the outer parts, as is observed in 
cD galaxies. 

Key words: galaxies: elliptical and lenticular, cD - dark matter - galaxies: interac- 
tions - galaxies: kinematics and dynamics - galaxies: structure 



1 INTRODUCTION 

Brightest cluster members (BCMs) are giant ellipticals 
found in the central parts of clusters of galaxies. Some of 
these galaxies are further classified as cD galaxies because 
they show in their external parts a projected surface bright- 
ness excess over the standard r^'** law. Several models have 
been proposed for the origin of BCMs in general and of cDs 
in particular. A first class of models relies on the presence 
of cooling flows in clusters of galaxies. If the central clus- 
ter density is sufficiently high, intra-cluster gas can gradu- 
ally condense and form stars at the bottom of the potential 
well (Cowie & Binney 1977; Fabian & Nulsen 1977). Such 
models, however, imply the formation of a lot of new stars, 
of which there is no good observational evidence. A second 
class of models relies on the dynamical mechanisms of merg- 
ing, cannibalism and accretion (Ostriker and Tremaine 1975; 



Ostriker and Hausman 1977; Hausman and Ostriker 1978; 
McGlynn and Ostriker 1980 etc.). Schematically the central 
giant galaxy is formed by the merger of smaller galaxies and 
its growth continues by further mergings and cannibalism. 
It also accretes material stripped from the outer parts of the 
other galaxies, by interactions between galaxies, or between 
a galaxy and the common background. It has been proposed 
(Gallagher & Ostriker 1972; Richstone 1976) that such ma- 
terial could be responsible for the luminous envelope of cD 
galaxies. Merritt (1984) argued that most of the formation 
happens during cluster collapse. 

Most of the mass of clusters or groups of galaxies is not 
bound to the galaxies that constitute them. Instead it is dis- 
tributed in a common background encompassing the whole 
cluster and it includes the high temperature gas emitting X- 
rays and an unknown fraction of dark matter. This unknown 
mass distribution must have influenced the formation mech- 
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anism of BCMs. Nevertheless, considerable uncertainty still 
remains as to what fraction of the total cluster mass resides 
in the dark matter, and how it is distributed. 

Observations of galaxy clusters can be used to constrain 
their mass distribution. The first attempts to obtain the 
profile of this mass distribution were based on galaxy counts 
and the projected velocity dispersion of the galaxy system. 
They rely on the assumption that the cluster has reached 
virial equilibrium and depend on what is assumed about 
the possible anisotropy of the system. They lead to mass 
distributions possessing central cores where the mass density 
is nearly constant (e.g. Kent & Gunn 1982). More recent 
studies, however, point to a mass distribution with a cusp 
in the central parts (Carlberg et al. 1997). 

For some galaxy clusters which are very luminous in X- 
rays, many authors have studied the total mass distribution 
under the assumption that the hot gas is in hydrostatic equi- 
librium within the total potential. The correlation between 
the velocity dispersion of the galaxies and the X-ray temper- 
ature (e.g. Wu et al. 1998), as well as the observational evi- 
dence that the dynamical properties of clusters have evolved 
little since z ~ 0.8 (Mushotzky & Scharf 1997; Bahcall, Fan 
& Chen 1997; Henry 1997; Rosati et al. 1998; Vikhlinin et al. 
1998, etc.) indicate that clusters are globally relaxed struc- 
tures. Further supportive evidence comes from the fact that 
Crimele et al. (1997) find a relation between the gas and 
galaxy densities consistent with the assumption of hydrody- 
namical equilibrium. Thus the Jeans equation together with 
the hypothesis of hydrostatic equilibrium can be used to de- 
scribe clusters globally, except for clusters which are still in 
the process of formation, or of merging with another cluster. 
In most cases X-ray mass estimates have been obtained with 
models with flat cores. More recently, however, Tamura et 
al. (2000) have shown that models of the total mass distri- 
bution with cusps give better flts to the data of the cluster 
A1060 than models with flat cores. 

Yet a third way to derive the mass in a cluster is to 
use the gravitational lensing effect which is produced when 
the images of background galaxies are distorted by the mass 
concentration of a galaxy cluster. In the case of strong lens- 
ing background galaxies are distorted as giant blue arcs. A 
massive foreground cluster can also change the shapes (shear 
effect) and the number density (magniflcation) of the faint 
background population. These effects are known as weak 
lensing. Arcs have been used intensively to map the matter 
distribution in the central parts of clusters (Fort & Mel- 
lier 1994; Mellier 1999 and references therein). In particular 
considerable progress was achieved from recent spectacular 
images of arc(let)s obtained with the HST. Many different 
mass distributions, including profiles with cores or cusps, 
have been used in the modelling. There seems to be agree- 
ment that, if cores exist, their radius should be less than 
50/ij~oo Kpc (e.g. Mellier 1999 and references therein). Ex- 
cept for that, no clear consensus seem to have been reached 
on the form of the radial profile in the innermost parts. 

Comparison of the results obtained with the various 
methods described so far lead to considerable controversy. 
In particular, results from strong lensing are a factor 2 — 4 
higher that results from X-rays. (e.g. Loeb & Mao 1994, 
Miralda-Escude & Babul 1995, Tyson & Fischer 1995; Wu 
et al. 1998, Melher 1999). This could be explained by the 
fact that, although the X-ray gas should be globally in grav- 



itational equilibrium, locally in the center-most parts virial 
equilibrium may not be achieved. It should, however, be 
noted that models with smaller core radii, with cusps, or 
more generally with higher masses in the center-most parts, 
reduce the differences between the results obtained with the 
different methods (e.g. Wu et al. 1998, Makino & Asano 
1999, Markevitch et al. 1999). 

N-body numerical simulations of hierarchical gravita- 
tional collapse of an ensemble of cold coUisionless particles 
lead to halo radial profiles with a cusp in the central parts. 
Navarro, Frenk and White (1996) argued that such profiles 
can be fitted by a "universal" profile of the form: 



p(r-) 



Ps 



(r/rs)(l + r/rs 



(1) 



where ps and r^ axe the characteristic density and length 
respectively. This profile was first obtained in the context 
of the standard cold dark matter cosmology, but subsequent 
numerical studies have shown that this profile is independent 
of halo mass, of the initial density fluctuation and of cos- 
mology (e.g. Cole & Lacey 1996; Navarro Frenk and White 
1997; Eke, Navarro & Frenk 1998; Jing 1999). Simulations 
with higher numerical resolution result in profiles which are 
steeper in the central parts, i.e. of the form p ~ r"^* — r~^'^ 
(Fukushige & Makino 1997, Moore et al. 1998), while cos- 
mologies using collisional dark matter point to even steeper 
(p ~ r"^) profiles (Yoshida et al. 2000; Moore et al. 2000), 
at least in the limit of short mean free paths. 

For these considerations it seems appropriate to com- 
pare the evolution of clusters of galaxies within backgrounds 
with different mass profiles, in hope that this may lead to 
clues about which kind of mass distribution is more realis- 
tic. Several N-body simulations of the dynamical evolution 
of large groups or clusters of galaxies with each galaxy being 
represented by many particles have been reported in the lit- 
erature. Funato, Makino & Ebisuzaki (1993), Sensui, Funato 
& Makino (1999) and Garijo, Athanassoula & Garcia-Gomez 
(1997, hereafter GAG) use model clusters where initially all 
the mass is in the galaxies. Funato, Makino & Ebisuzaki 
(1993) argue that the Faber- Jackson relation is a result of 
the evolution of galaxies driven by interaction with other 
galaxies and with the tidal field of the parent cluster. Sen- 
sui, Funato & Makino (1999) follow the evolution of a cluster 
of 128 galaxies initially in virial equilibrium with all mass 
being attached to the galaxies. They find that within a few 
cluster crossing times half of the total mass has escaped indi- 
vidual galaxies, and that the amount of stripped material is 
larger for galaxies in the central regions. The density profile 
of the common halo thus formed is cuspy, with a radial de- 
pendence roughly p ~ r~^ in the central region and p ^^ r'"^ 
in the outer regions. Bode et al. (1994), GAG and Dubin- 
ski (1998) witnessed in their simulations the formation of 
a massive central galaxy. Bode et al. (1994) modelled poor 
clusters with 50 galaxies and a varying percentage of the 
mass in a common background. They find multiple nuclei 
in the central object in between 10 and 40 per cent of the 
time. Increasing the percentage of mass in the common back- 
ground slows the merging rate and this can be sufficient for 
stalling the formation of the central object if 90 per cent 
of the mass is in a common background. GAG studied in 
detail the properties of the central object and found good 
agreement with those of brightest cluster members. The ob- 
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ject formed is roughly spherical, oblate or mildly triaxial. In 
the latter cases the orientation of the central object corre- 
lates well with that of the central group. The triaxiality is 
in general stronger in the outer parts of the central object. 
They also find three types of projected density distributions. 
Objects in the first category have profiles well fitted by an 
r^'^ law. In the second category of objects the radial profile 
follows the r^'* law only in the main part of the galaxy and 
falls bellow it in the outer parts. In the third, most inter- 
esting, category the projected surface density of the objects 
follows the r^'*' law in the main body of the galaxy, but 
is systematically above it in the outer parts, as is the case 
for cD galaxies. The projected velocity dispersion profiles of 
the central objects also agree with those of brightest cluster 
galaxies. Dubinski (1998) used an N-body simulation of a 
cluster of galaxies in a hierarchical cosmological model. He 
replaced 100 dark matter halos at 2 = 2 with self-consistent 
disc-|-bulge-|-halo galaxy models, thus attaining high resolu- 
tion. He witnessed formation of the central galaxy through 
the merging of several galaxies along a filament. The central 
object was flattened and triaxial, showing alignment with 
the primordial filament. Its projected surface density was 
well described by an r^'* law and showed no excess in the 
outer parts, as one would expect from a cD. The kinematics 
of the object also was in good agreement with the observa- 
tions. 

In this paper we will study the influence of the dark 
matter distribution on the evolution of the cluster and on 
the formation and properties of the central galaxy. Our aim 
is to show that the observable properties of this galaxy can 
be used to constrain the dark matter distribution in the cen- 
tral parts. The paper is organized as follows. In section 2 we 
present the simulations and the initial conditions. In sec- 
tion 3 we describe the different evolution of the various 
models and in section 4 we show how the different initial 
conditions influence the formation process and the proper- 
ties of the central giant galaxy. Section 5 discusses the main 
assumptions and simpliflcations entering in our simulations. 
Our results are discussed and summarised in Section 6. 



2 SIMULATIONS 

We have performed a total of 17 simulations following the 
time evolution of a group of 50 galaxies in virial equilibrium 
within a common halo. In nearly all cases we used a total 
of 100 000 particles with the same mass. A fraction of these 
particles was used to represent the galaxies and the rest were 
used for the common halo. In one simulation only one third 
of the total mass was in galaxies and for this case we used 
200 000 particles, while maintaining the total mass constant, 
in order to have galaxies represented by a reasonable number 
of particles. 

The common halo was represented alternatively as a 
Plummer or as a Hernquist sphere (Hernquist 1990). The 
density distribution of a Plummer sphere follows the law: 



Table 1. Initial conditions 



ppix) = 



3Mc,, 
47raS 



+ fflp 



5/2 



(2) 



Plummer models 
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CLp 
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Npa 


Common Halo 


P5f3/2rl 
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3/2 
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1000 


Live 


P15f3/2rl 
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P10f3/2rl 


10 
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1000 
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PlOflrl 


10 


1 
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P10f3/4rl 


10 
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1000 


Rigid 


P10flr2 


10 
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1333 


Live 


P10flr4 


10 




4 


1600 


Live 


P10flr6 


10 
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1714 


Live 


PlOflrlO 


10 




10 


1818 


Live 


P10flr20 


10 




20 


1904 


Live 


PlOrl 


10 


00 
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1000 
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P10r20 


10 


00 


20 


1904 


Live 


Hernquist models 


Name 


ah 


/ 
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Nva 


Common Halo 


HlOflrl/2 


5.65 
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0.5 


1333 


Live 


HlOflrl 


5.65 


1 
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1000 


Live 


H10flr5 


5.65 
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1666 
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H10flr20 


5.65 
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20 


1904 


Live 



which represents the scale-length of the distribution. The 
Hernquist distribution follows the law: 

Mch ah 1 i„s 

Ph{r) = — -f-— Tg. (3) 

2-K r [r + ahy 

Two parameters are also necessary to determine this dis- 
tribution: the total mass in the common halo A'IcH and the 
scale- length a^- Note that while the Plummer distribution 
has a central density of SMch/^Tvap, the Hernquist distribu- 
tion is singular at the origin. Its projected surface density 
follows closely the r^''^ law for a large fraction of the radius 
(Hernquist 1990). It gives a cusped density distribution in 
the central parts which falls as p{r) ^ r~^ for small radii 
and as p{r) ~ r^* at large radii. 

For a given mass of the common halo Mch, the param- 
eters Up and ah were selected as to give a common halo 
with the same total energy. This gives the following relation 
between the two parameters: 



ah 



16 
— c 
97r 



(4) 



This distribution is characterised by the parameters Mch, 
which represents the total mass in the common halo, and Op 



In this way we can compare simulations with the same global 
properties changing only the mass profile of the common 
halo. 

The positions and velocities of the particles in the com- 
mon halo are drawn according to the chosen Plummer (cf. 
eg. Aarseth, Henon & Wielen 1974) or Hernquist model 
(Hernquist 1990) . The positions and velocities of the center 
of masses of the galaxies are drawn from the same distribu- 
tion as the common halo, except that we impose a cut-off 
radius which is a fraction of the half mass radius R1/2 of 
the common halo. Thus the galaxies are distributed within 
a radius selected as 



Rci = fRi/2 



(5) 
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where / is a parameter which is used to determine the initial 
concentration of the galaxy system. Finally, we also vary 
the ratio between the mass in galaxies and the mass in the 
common halo. This mass ratio is given by the parameter: 






(6) 



where Ng represents the initial number of galaxies and Mg 



the initial mass of each galaxy. We use in all cases Ng 



50 



identical galaxies, whose density profiles follow a Plummer 
law with a scale-length of Ugai = 0.2. The mass of each 
galaxy, Mg, will depend on the mass ratio selected for each 
particular initial condition, i.e. on the value of r. 

We repeated one of the simulations using a rigid halo 
which is not modified during the whole simulation, instead 
of a live halo composed of particles. We did this in order 
to study the effect of dynamical friction and the possible 
differences introduced by a halo capable of responding to 
the variation of the mass distribution as the central giant 
galaxy is formed. 

In Table 1 we give more detailed information about our 
simulations. In Column 1 we list the simulation name, coded 
conveniently so as to describe the initial conditions of each 
simulation. The first letter P or H is used to note the type of 
mass distribution, P for Plummer models and H for Hern- 
quist models. This letter is followed by one or two digits 
indicating the scale-length of the common halo in simula- 
tion units. For the case of Hernquist models we indicate 
the scale-length of the Plummer model with the same to- 
tal energy. Thus from the name of the simulations we can 
readily see which of the different simulations are directly 
comparable. After these digits comes the letter / followed 
by the value of this parameter and finally the letter r fol- 
lowed by the value of the mass ratio. In this work we have 
considered values of r between 1/2 and 20, thus covering a 
range of more than one order of magnitude. The three val- 
ues of the parameters determining the initial conditions are 
repeated for clarity in Columns 2, 3 and 4. In Column 2 for 
the case of the Hernquist simulations we list the value of the 
scale-length of the Hernquist distribution. In Column 5 we 
indicate the total number of particles used to represent each 
galaxy and finally in Column 6 we indicate if the common 
halo is rigid or a live one evolving freely. In two of the simula- 
tions, PlOrl and PlQr2Q, the initial positions of the galaxies 
were selected anywhere within the common mass distribu- 
tion without any cut-off. Finally in the Hernquist simula- 
tion with low galaxy-to-common-halo mass ratio, namely 
HlOflrl/2, the total number of particles was doubled, while 
the total mass was kept constant. 

To evolve these systems we used the GRAPE-3AF 
system of the Marseille Observatory using a version of a 
treecode specially adapted for this system (Athanassoula et 
al. 1998). This version, since the tree is descended for a 
whole group of particles, rather than for every particle sep- 
arately, is more precise than a standard treecode (Barnes & 
Hut 1986). For the opening angle we used 9 = 0.7, for the 
time step 5t — 0.0078125 and for the softening parameter 
e = 0.0625. With these values the total energy was con- 
served to within a few parts in a thousand. The simulations 
with a Plummer matter distribution were followed for a to- 
tal of roughly 60 computer units or, equivalently, roughly 
7680 time steps. Each simulation required about 24 hours of 



CPU time on the Marseille five board GRAPE-3AF system. 
The evolution of simulations with a Hernquist distribution 
was followed for twice as long, to follow best the evolution 
of the central object. 

In this paper, we use the same units as in GAG, where 
the simulation units were selected as G = 1, the total mass 
Mt = 50 and the typical scale-length of a galaxy Ogai = 0.2. 
With the aim of comparing with the observations we set 
Mt = 1.5 X 10^'^Mq and Ogai ~ 6 Kpc. These values give 
a total of 1.4 x 10 years for the time unit, 210 Km/s for 
the velocity unit and a total simulation time of T = 8.4 and 
16.8 Gyr, for Plummer and Hernquist distributions respec- 
tively. Note that this selection of conversion factors, albeit 
reasonable, is not unique and neighbouring values could also 
be used. This should be taken into account when comparing 
our simulations with observational data, hence agreements 
to within a factor of two should be considered to be quite 
satisfactory. 



3 EVOLUTION OF THE SIMULATIONS 

3.1 A general impression from snapshots 

As in the simulations without common halo (GAG) the most 
noticeable result of the dynamical evolution is the formation 
of a giant central galaxy in the central parts of the group. 
The rate of formation of this galaxy, as well as its final ex- 
tension, depend on the properties of the common halo and of 
the galaxy system. This can be seen qualitatively in Figs. ^ 
to 0, where we show respectively the influence of the scale- 
length, of the initial extension of the galaxy distribution, of 
the ratio of the mass in galaxies to the mass in the common 
halo and the influence of different mass profiles. 

Amongst our simulations with Plummer background 
the one that starts more centrally concentrated is P5f3/2rl 
and the one that starts off with the least central concentra- 
tion is P15f3/2rl. Since the other parameters are the same 
in the two simulations, the comparison of the two evolu- 
tions, given in Fig. hi, shows the effect of the scale-length 
of the galaxy distribution and of the common halo, which 
changes by a factor of three between the two cases. The more 
concentrated systems, with shorter scale-length, have faster 
evolution rates. In this way in simulation P5f3/2rl a central 
galaxy is rapidly formed which controls thereafter the dy- 
namical evolution of the system, merging with the satellite 
galaxies in a shorter time than in the rest of the simulations. 
On the other hand, in the simulation P15f3/2rl no central 
giant galaxy is formed and practically all the galaxies sur- 
vive until the end of the simulation. 

In Fig. Q we compare the evolution of simulations 
P10f3/2rl and P10f3/4rl to understand the effect of chang- 
ing the initial cut-off radius of the galaxy distribution. Since 
the number of galaxies is initially the same in both cases 
a decrease of the cut-off radius entails a higher density of 
the galaxy system. As we can see, if the galaxy distribu- 
tion is initially denser, the system has a faster evolution 
rate. Since this is also the effect of a more concentrated 
distribution, as discussed earlier, we can conclude that it 
is the density in the central parts that controls the rate of 
evolution. Thus distributions with higher central densities, 
coming either from higher central concentrations or higher 
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Figure 1. Comparison of the evolution of two simulations, one corresponding to a small scale-length of the common halo and galaxy 
distribution (left panels) and the other to a large one (right panels). Three times are shown for each simulation, the initial one, one not 
far from half way through the simulation and the third one at the end. The times are given in the upper right corner of each panel and 
the simulation name in the top panel of each column. Only one particle in 20 is plotted. 
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Figure 2. Effect of tfie initial cut-oflf radius of the galaxy distribution on the evolution of the group. If the cut-off radius is smaller, i.e. 
if the galaxy system is initially more dense (right panels) then the evolution of the system is faster. The layout of this figure follows that 
of Fig. §. 
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Figure 3. Effect of the ratio of tfie mass in galaxies to the mass in the common halo on the evolution of the group. If the mass is initially 
placed mainly in the galaxies (right panels) the system evolves faster than if a larger fraction of the mass is in the common halo (left 
panels) . The layout of the figure follows that of Fig. hi 
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Figure 4. Effect of the diflfcrent common halo and galaxy mass distribution. The left panels show the evolution of the galaxy system 
within a Hernquist model while the right panels shown the evolution of the galaxies within a Plummer model of the same total energy. 
The central galaxy grows faster in the Hernquist models. The layout of the figure follows that of Fig. hi 
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densities over a larger area, will have a faster evolution than 
distributions with lower densities in the central area. 

Fig. H compares the evolution of simulations PlOflrl 
and P10flr20 and shows the effect of the ratio of the mass 
in individual galaxies to the mass in the common halo. Simu- 
lation PlOflrl has 20 times more mass in the common halo 
than simulation P10flr20. The comparison shows clearly 
that the simulation with the lower fraction of the mass in 
galaxies (PlOflrl), i.e. with the higher fraction of the mass in 
the common halo, evolves slower than the simulation with 
the higher fraction of the mass in galaxies. This in agree- 
ment with what was seen in the simulations of Bode et al. 
(1994) for poor clusters of galaxies, as well as those of Barnes 
(1985), Bode et al. (1993) and Athanassoula et al. (1997) for 
compact groups. This may be easily understood (Athanas- 
soula et al. 1997) if we think of the galaxies as a perturbation 
of the global common halo potential. If the fraction of the 
mass in galaxies is higher, the perturbation is also stronger, 
the galaxies will attract each other more strongly and have 
more chances of colliding and therefore of merging. The op- 
posite will be true if the fraction of the mass in galaxies is 
small. In the limiting case were the galaxies are considered 
as test particles in the common halo, they will collide only 
if their orbits accidentally cross. 

Fig. W compares the evolution of two simulations, one, 
HlOflrl, in which the distribution of the common halo and 
the galaxies follows a Hernquist law, and the other, PlOflrl, 
in which it follows a Plummer law, the remaining parameters 
being the same. The difference is quite striking. We note that 
the evolution is much more rapid for the Hernquist model. 
Thus, at any given time, for the Plummer model there are 
more galaxies that have not fallen into the central one than 
for the Hernquist model. 

3.2 Global parameters 

We will use in the analysis of the simulations the same def- 
initions of the central galaxy and of the merging conditions 
as in GAG. They allow us to make a more quantitative study 
of the evolution of the system. The first factor that we can 
consider is the number of galaxies that survive at a given 
time. This is plotted as a function of time in Figures H and 
In the upper panel of Fig. we show the results for the 
simulations with a Hernquist distribution and in the lower 
panel for the simulations with a Plummer distribution. In 
this figure we concentrate on simulations with Op = 10 and 
/ = 1 and vary only the fraction of mass in the galaxies. 
The evolution shown in Fig. f\ is in good agreement with the 
qualitative results of Fig. pi In general a larger percentage of 
the total mass in galaxies leads to a faster evolution, for the 
reasons explained in the previous section. For simulations 
with Plummer distributions the effect is only important in 
the second half of the evolution when the central galaxy is 
sufficiently massive. On the other hand for simulations with 
Hernquist distributions the effect is clear from the start, and 
is stronger. These differences are due to the different evolu- 
tionary histories of the two types of simulations, which will 
be discussed in section 4.1. Finally we should also note that 



i 





PlOflrSO 
PlOflrlO 

. P10flr4 
. PlOflrS 
. PlOflrl 



the fastest of all evolutions is for simulations with Hernquist 
profiles and the biggest fraction of the mass in galaxies. 

Fig. B also shows the number of galaxies that have not 
merged at a given time as a function of this time, but concen- 



Time (Gyr) 

Figure 5. Time evolution of the number of surviving satellite 
galaxies. Comparison of the effect of the initial fraction of mass 
in galaxies both for Hernquist (upper panel) and Plummer (lower 
panel) mass distributions with the same total energy. 



trates on Plummer profiles, to see the effect of the three free 
parameters of the Plummer distribution. The upper panel 
illustrates the effect of the scale-length of the distribution 
of the common halo and of the galaxies, the middle one the 
effect of the cut-off radius and the lower one the effect of 
the fraction of mass in the galaxies. The results agree well 
with what is seen qualitatively in Figures hi to M. The galaxy 
number decreases faster in the cases of more concentrated 
mass distribution, i.e. halos and galaxy distributions with 
a smaller scale-length. Also configurations with a smaller 
cut-off radius of the galaxy distribution evolve faster. 

The next quantity that we can consider is the evolution 
of the mean radius of the system of surviving galaxies. It 
is calculated as the average distance of the centers of the 
surviving galaxies from the center of the group not taking 
into account the giant central galaxy. Its time evolution is 
shown in Fig. fTl As in Fig. fa we show in this figure the sim- 
ulations with ap = 10 for the Hernquist halo distributions 
(upper panel) and Plummer distributions (lower panel). In 
all but a couple of cases the mean radius of the system stays 
constant. This argues that the satellite galaxies that fall into 
the central one are not only those that have lower binding 
energies and orbits with a smaller mean radius. Rather all 
galaxies, independent of their binding energy, are concerned. 
The case of H10flr20 is different, since the mean radius of 
the surviving galaxies increases with time. This is due to the 
fact that one galaxy is on an outwards going orbit, while 
the total number of galaxies is rapidly reduced (cf. Fig. |5|) 
so that this one galaxy can inffuence the statistics. 

The mean mass in surviving galaxies, i.e. not taking into 
account the giant central galaxy, is indicative of whether the 
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Figure 6. Time evolution of the number of surviving satellite 
galaxies. Comparison of the effect of the different parameters 
of the initial configuration. The upper panel shows the effect of 
the scale-length, the middle panel the effect of the cut-off in the 
galaxy distribution and the lower panel the effect of the fraction 
of mass in galaxies. 
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Figure 7. Time evolution of the mean radius of the galaxy svstem 
of surviving galaxies. The layout of the figure is as for Fig y. 
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Figure 8. Time evolution of the mean mass of the survivi: 
galaxies in computer units. The layout of the figure is as for Fig 



mass loss effects due to stripping are important in the dy- 
namical evolution of the groups, or not. The time evolution 
of this parameter is shown in Fig. H for Hernquist systems 
(upper panel) and for Plummer systems (lower panel). In all 
the cases the mass in galaxies diminishes with time, which 
means that stripping is present and that the stripped mass 
is not all captured by other individual galaxies. Some of it 
is captured by the giant central galaxy, while another part 
stays in between the individual galaxies. 



3.3 The rigid halo case 

If we could, in general, use in simulations a rigid halo, ren- 
dered by an external forcing, constant in time, instead of a 
"live" halo formed by a large number of particles, we could 
save a considerable amount of CPU time. For this reason 
we will, in this section, compare the effects of a rigid and a 
live halo on the evolution of a simulation, to see whether the 
corresponding gain in CPU is possible without introducing a 
bias in the results. At the beginning of the simulation a live 
halo will give, on average, the same forces as a rigid halo. 
This will, however, not be necessarily the case as the simu- 
lation evolves, since the mass distribution of the halo may 
change with time, in response to the evolution of the dis- 
tribution of the individual galaxies and to the formation of 
the giant central galaxy. Furthermore, even in cases where 
the mass distribution in the common halo hardly evolves, 
we are not assured of the adequacy of the rigid halo, since 
it will not reproduce the dynamical friction on the satellite 
galaxies that the live halo exerts naturally. This effect could 
in principle be estimated by the Chandrasekhar (1943) for- 
mula, except for the fact that the common halo does not 
necessarily have a homogeneous density distribution, while 
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Figure 9. Effect of the rigid halo on the evolution of the global 
properties of the system. The upper panel compares the number 
of surviving galaxies, the middle one their mean radius and the 
lower one their velocity dispersion. 
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Figure 10. Effect of a cutoff in the initial distribution of the 
galaxies on the global properties of the system. The upper panel 
compares the number of surviving galaxies and the lower one their 
velocity dispersions. 



dynamical friction is much more important in strongly cen- 
trally concentrated backgrounds than in less concentrated 
ones (Athanassoula, Makino & Bosma 1997). It is thus nec- 
essary to compare a simulation with a live halo to one with 
a rigid one, to see how big a bias the rigidity of the halo 
introduces. 

To check this effect we have repeated one simulation, 
namely P10f3/4rl, using a rigid halo instead of the live halo 
composed of 50000 particles. In Fig. we compare the evolu- 
tion of the global parameters of the two systems. We can 
see that the number of surviving galaxies is considerably 
lower in the case of a live halo (upper panel). The principal 
reason for this is the fact that the dynamical friction with 
the live halo brakes the galaxies and, as they move slower, 
they merge easier with the central galaxy. The radius which 
contains half the mass of the galaxy system stays in both 
cases roughly constant with time and the average value is 
roughly the same in the two simulations (middle panel) . The 
velocity dispersion of the system of galaxies shows in both 
cases an increase with time, presumably due to the fact that 
the faster individual galaxies merge less than the slower one. 
The lower panel of Fig. H shows that this increase is stronger 
for the live halo than for the rigid one. This should be linked 
to the fact that there are more mergings in the live halo case. 



3.4 The effect of the cutoff on the galaxy 
distribution. 

In most of our simulations we have included an artificial 
cutoff in the initial distribution of the galaxies. This partly 
reflects the fact that the distribution of galaxies in a group 



or cluster also does not exceed to infinity. It has the added 
advantage of improving the statistics for a fixed number of 
particles per simulation. In this section we discuss its effect 
by comparing the evolution of runs with cut-off with that 
of runs without cut-off and otherwise the same properties. 
Fig. nG compares two such sets of simulations, simulation 
P10r20 to P10flr20, and simulation PlOrl to PlOflrl. In 
this figure we compare two global properties of the system 
of surviving galaxies: their number (upper panel) and their 
velocity dispersion (lower panel). The number of surviving 
galaxies decreases considerably faster with time if we include 
a cut-off. This is due to the fact that by constraining the vol- 
ume, while keeping the number of galaxies fixed, we increase 
the density of the galaxies and this in turn, as discussed in 
the beginning of this section, leads to a faster evolution of 
the system. 

The lower panel shows that the velocity dispersion is 
smaller for the cases without cut-off. There can be two pos- 
sible explanations to this. The first one relies on the fact that 
it is the galaxies with the lower velocities that will merge, 
and therefore systems that involve more merging should also 
have a higher dispersion of velocities. If, however, this was 
the explanation then we would note differences in the veloc- 
ity dispersions after the differences in the number of remain- 
ing galaxies occurred, which, as can be seen in Fig. nil is not 
the case. We will thus explain the difference by the fact that 
the velocity dispersion of the group decreases with radius. 
Thus a more extended group has a lower mean dispersion of 
velocities than a more compact one. 
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Figure 11. Time evolution of tlie total mass in the central parts, 
excluding the mass in individual galaxies. Left panels correspond 
to Hernquist simulations and right panels to Plummer simula- 
tions. The solid line corresponds to the evolution of the total 
mass while the dotted line gives the fraction of this mass that 
comes from tidal stripping. The name of each simulation is given 
in the upper left corner of each panel. 



4 THE CENTRAL OBJECT 

The central object is defined as in GAG. The properties of 
this object are studied using the same methods described in 
that paper and we refer the reader to that paper for more 
details. Here we will discuss directly the physical properties 
of the objects formed in our simulations. 



4.1 Stripping versus merging 

The mass of the central galaxy grows with time by means 
of two mechanisms: by merging of satellite galaxies and by 
accreting material tidally stripped from galaxies. We are in- 
terested in how the relative importance of these two mecha- 
nisms depends on the initial conditions. In Fig. [nl we show 
for some selected simulations the time evolution of the total 
mass in the central regions not belonging to the individual 
galaxies. In the left panels we show the simulations where 
the initial density distribution follows a Hernquist profile, 
while in the right panels we show the cases of initial Plum- 
mer density distributions. Each horizontal couple of panels 
corresponds to two simulations with identical or compara- 
ble global properties. The solid line corresponds to the total 
mass as a function of time, while the dotted line corresponds 
to the part of this mass that comes from stripped material. 
The difference between the two lines at a given time gives 
the part of the total mass that comes from mergers with 
satellite galaxies. 

Let us first note that in general the total mass of the 



central object is larger in simulations starting with Hern- 
quist profiles than in simulations starting with Plummer 
profiles. The total mass of the central object (at a given 
time) is also larger in simulations where the total mass that 
is initially in galaxies is larger, as expected. 

As can readily be seen from Fig. hll the initial distribu- 
tion of the mass profile (background mass and distribution 
of galaxies) is determinant not only for the total mass of the 
central object, but also for the formation process leading to 
it. In Plummer simulations all the initial mass increase is 
due to the accretion of stripped material. This is true for 
a longer period of time in the case of simulation PlOflrl, 
where the mass in galaxies was initially the same as the mass 
in the common background. In this simulation also the final 
fraction of the total mass that comes from tidal stripping is 
more important than in the others. On the other hand in the 
Hernquist simulations the merging of satellite galaxies into 
the central galaxy is present from the start. This should be 
due to the differences in central concentration. Our results 
argue that the central concentration of the Hernquist pro- 
file is sufficiently high to compensate for the higher relative 
velocities of the galaxies and to cause a high initial burst of 
mergings. 

Thus, as can be seen from Fig. hll the relative impor- 
tance of the stripping is more significant in the Plummer 
cases than in the Hernquist ones. Conversely the absolute 
amount of stripped material in the first part of the simula- 
tion is larger for the Hernquist cases. This difference is due 
to the fact that the total mass of the central object in the 
first part of the simulation is considerably bigger for Hern- 
quist cases. 

These differences between simulations with Plummer 
and with Hernquist profiles can be readily understood from 
the different evolutionary histories of the two groups of simu- 
lations. In the Plummer case the mass distribution in the 
central parts is nearly flat and as the initial position of the 
satellite galaxies is selected following this distribution, no 
initial concentration of galaxies is present. Thus, there is no 
strong merging of satellite galaxies into a central object dur- 
ing the first steps of the simulations. The satellite galaxies 
orbiting within the core radius occasionally merge with some 
of their companions. These mergings, which do not neces- 
sarily occur in the center-most parts increase the number of 
larger galaxies within the core radius, until two or more of 
these finally merge to give a giant central object. 

On the other hand, the dynamical evolution of the simu- 
lations with a Hernquist profile is very different. The strong 
peak of this density profile towards the central parts gives a 
strong initial concentration of satellite galaxies which merge 
very fast to create the seed of a giant central object. We 
refer to this stage of the evolution as an initial burst of 
mergers. The central object will accrete some stripped ma- 
terial and merge with some other satellite galaxies entering 
the central parts. After the strong initial merger, the mass 
increase of the galaxy is milder than in the initial steps of 
the simulations. As we will see in subsection 1.3 these two 
different formation stories give galaxies with very different 
observational properties. Not all the central objects will have 
profiles that resemble real giant elliptical or cD galaxies. 
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Figure 12. Volume density profile, given by filled circles, as a 
function of radius for simulations PlOflrl (upper left), P10flr20 
(upper right), HlOflrl (lower left) and H10flr20 (lower right). 
The time is given in the upper right corner of each frame. The 
solid lines give eyeball fits by straight lines. 



4.2 Volume density profiles. 

Figure |l2| shows the volume density profile as a function 
of radius for simulations PlOflrl, P10flr20, HlOflrl and 
H10flr20 at t = 8.2 Gyrs. Plummer simulations show clearly 
two distinct regions, the outer one being much steeper than 
the inner one. Each of these can be relatively well fitted by a 
straight line in a log- log plane. The slopes of these two lines, 
as obtained by eyeball fits, are roughly -2 and -5. The outer 
part is clearly discernible all through the simulation and its 
slope, which is roughly equal to the slope of the outer part 
of the Plummer density profile, stays roughly constant all 
through the simulation. This is not the case for the inner 
part, which in some earlier times is more difficult to discern 
than at i = 8.2 Gyrs. Since our simulations stop at T = 8.4 
Gyrs, we have no information on the evolution of the profile 
at later times. 

The density profiles coming from Hernquist simulations 
do not show the clear existence of two distinct regions, al- 
though they do steepen at larger radii. For consistency we 
also made eyeball fits with two straight lines. The inner 
slopes thus found are steeper than the corresponding ones 
for the Plummer simulations, while the slopes of the parts 
further out are roughly the same. The profile of H10flr20 
has a shallower sloper in the outermost parts, which will be 
discussed at length in the next subsection. 

We compared all these profiles to the Hernquist profile 
(Hernquist 1990) and found that, if we exclude the inner- 
most and outermost parts, the Plummer simulations with 
little mciss in the common halo, as well as all Hernquist 
simulations, are well fitted by this profile. This was, how- 



ever, not the case for the Plummer simulations with a high 
fraction of the mass in the common halo. The physical im- 
plications of this result will be discussed further in the next 
subsection with the help of the projected density profiles. 

4.3 Projected density profiles. 

The most interesting property of the central galaxies formed 
in these simulations appears when we study their projected 
density distributions. Remnants of the merging of a pair of 
progenitor galaxies have a projected density profile which 
is well fitted, in all but the innermost parts, by the r^'* 
law (see for instance Barnes and Hernquist 1992, or Barnes 
1998 and references therein). A similar result was found for 
multiple mergings in compact groups (Weil & Hernquist 
1996, Athanassoula & Vozikis 1999). Furthermore the ob- 
jects formed in GAG from multiple mergers in large groups 
or small clusters also showed this property, except for an 
excess projected density in the outermost parts in some cD- 
like cases. Our new simulations with a common halo present 
a larger variety of profiles. 

The three-dimensional density profiles of the central ob- 
ject formed in the Plummer models with a considerable frac- 
tion of the mass in the background halo can not be well fitted 
by a Hernquist law, and therefore does not give an r^''^ in 
projection. In fact the larger the fraction of mass in the com- 
mon halo, the worst the agreement is. On the contrary, if the 
mass in the common halo is very much smaller than the mass 
in common galaxies (e.g. in simulation P10flr20), then we 
recover initial conditions closer to those of the simulations 
in GAG, and the profiles are better fitted by a Hernquist 
law. This form of the profiles does not depend either on the 
initial scale- length of the mass distribution, or on the ra- 
dius initially containing all galaxies, and depends only on 
the initial mass fraction in galaxies. These results are best 
illustrated in Figures flS^ and h4. 

The mass distribution in the central parts of simula- 
tion PlOflrl, where matter is equally distributed between 
the galaxies and the background, is heavily influenced by 
the matter in the common halo, and the violent relaxation 
which would have otherwise led to an r^'* projected pro- 
file is hampered. This result can be seen in Fig. [13, where 
we show the projected density profile of the central object 
formed in simulation PlOflrl at four different times. In each 
of the panels we show the result of 10 projections from ran- 
dom viewing angles, thus making clear that this result is 
projection independent. At no time does this profile follow 
an r^''^ law and thus the central object cannot be associated 
to any giant elliptical or cD galaxy, as are observed in the 
central parts of galaxy clusters. 

The central object formed in simulation P10flr20 feels 
only a small influence from the common halo and develops 
a density distribution which is well fltted by an r^'** law. 
This can be seen in Fig. [14, where we show the projected 
density profile of the central object at four different steps of 
the simulation. This effect is again not due to a particular 
viewing angle. As was already mentioned, this result is very 
similar to what was found in GAG since in P10flr20 there 
is very little inffuence from the background. 

Simulations where the common background is described 
initially by a Hernquist distribution give very different re- 
sults. The projected density profile for the central galaxy in 
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Figure 13. Projected density profiles of the central object formed 
in simulation PlOflrl at four different times. It is clear that the 
profiles are not well described by r^' * laws. In each panel we give 
by solid lines the result of 10 random projections. The time is 
given in the upper right corner of each panel and the projected 
density is in arbitrary units. 
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Figure 14. Projected density profiles of the central object formed 
in simulation P10flr20 at four different times. It is clear that the 
profiles are well described by r^'"* laws. The layout is as for the 
previous figure. 



simulation HlOflrl is shown in Fig. hsl As in the case of 
simulations PlOflrl, only half the mass was initially in the 
galaxies. We note that the projected density profile of the 



central object follows an 



nl/4 



law in the main body of the 



object and can thus be associated to a real giant elliptical. 
Simulations HlOflrS and H10flr20 have a much smaller 
fraction of their mass in the common halo; roughly 17 and 5 
per cent respectively. Their projected radial density profiles 
are shown for six different time steps in Figures hq and [ITI 
They show clearly that the radial profiles follow well the 
r^'*' law over the main body of the object. Here it can not 
be claimed that this profile is formed as a response to the 
massive common halo with a Hernquist radial profile, since 
the mass in the halo is rather small. These central objects 
were formed in a situation where there is a strong central 
concentration of galaxies at the start, thus favouring a strong 
and rapid merger. Thus giant galaxies formed from a rapid 
multiple merger in the central part of groups or clusters 
should have an r^'^ projected density profile over their main 
body. 



Now let us turn to the excess mass over the r 



1/4 



law 



in the outermost parts of the central object. This should be 
associated with the corresponding excess hght in the outer- 
most parts of some BCM's and giant elliptical galaxies. This 
light excess has, unfortunately, been called "halo", a term 
which could bring confusion with the dark halo. To avoid 
such a confusion, while being coherent with the nomencla- 
ture used so far, we will use the term "luminous halo" . Giant 
ellipticals with such a luminous halo are called cD galax- 
ies. We will now examine which of our simulations present 
central objects with such a luminous halo and which can 
therefore be associated with cD galaxies. 

Simulation P10flr20, which is the only one of those with 
an initial Plummer mass distribution that shows an r^'* 
profile over the main body of the galaxy, has only a very 
small density excess over the r^'* profile in the outermost 
parts of the central object, which hardly grows with time. 
It can thus not be easily associated with a cD galaxy. 

Central objects formed in simulations with a Hernquist 
background present such an excess of matter over the r^'* 
profile in the outermost parts. This is small in cases with 
a large fraction of the mass in the common background, as 
HlOflrl, and becomes more important for simulations where 
the fraction of mass in the galaxies is also more important. 
Thus the central object in simulation H10flr5 has a clear 
outer mass excess and that in H10flr20 an even bigger one 
(cf Figures [Lw and n/n), independent of the projection se- 
lected. The existence of this extra mass is clear relatively 
early on in the simulations, and increases considerably with 
time. Thus such objects should be associated to cD galax- 
ies, since they form in the central part of the group, follow 
an r^'* projected density law over the main part part of 
their body and have excess matter over that profile in the 
outermost parts. 

This inverse correlation between the excess mass in the 
outermost parts and the amount of mass in the common 
background halo could be due to the fact that a deviation 



from the 



rl/4 



more difficult to form in the strong 



.1/4 



background. Alternatively it could be due to the fact that 
the rate of evolution is slower in simulations with a massive 
common halo, i.e. that at a given time the three simulations 
are at different stages of their evolution. The total evolution 
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Figure 15. Projected density profiles of tlie central galaxy 
formed in simulation HlOflrl at six different time steps. The pro- 
files are well described by r^' ^ laws in the main body of the object 
and show luminosity excesses over the r^' ** in the external parts. 
In each panel we give by solid lines result of 10 random projec- 
tions and by a dashed line we repeat one of the profiles of the 
upper left panel. The time is given in the upper right corner of 
each panel and the projected density is in arbitrary units. 



time in these simulations was twice as long as that of Plum- 
mer simulations in order to study the possible influence of 
the dynamical evolution of the group in the distribution of 
the luminous halo material of the central object. In Figs. PTsl 
flM and nj we show that indeed the luminous halo of the cen- 
tral object increases with time, specially in the cases where 
there is a considerable fraction of mass initially in the galax- 
ies. As we have already pointed out in section 2, the con- 
version from computer units to physical units is not unique 
and thus the times given in figures ha, hq and hj could be 
modified if e.g. a different value was adopted for the masses 
or the scale-length of the initial galaxies. We nevertheless 
believe that all simulations were evolved sufficiently long for 
a clear picture to be obtained. 



5 LIMITATIONS OF OUR MODELS 

Simulations of the evolution of galaxy clusters, seen the com- 
plexity of the objects to be modelled, use a number of sim- 
plifications and assumptions, and the simulations presented 
here are no exception. We will here discuss the simplifica- 
tions and assumptions specific to our simulations to get a 
better understanding of how they might have influenced our 
results. 

The initial number of galaxies in all our simulations 
is 50. Thus they do not apply to rich clusters, but more 



Figure 16. Same as the previous figure but for simulation 
HlOflrS. 
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Figure 17. Same as in figure 15 but for simulation H10flr20. 
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to large groups, or poor clusters. They can also apply to 
sub-condensations in rich clusters - provided of course the 
remaining part of the cluster does not significantly influence 
the evolution - or subunits which will come together to form 
a large cluster. 

Our initial conditions are not cosmological. They were 
constructed specifically so as to allow us to study best the 
physics of the formation of the central object and partic- 
ularly the effect of different physical parameters on this 
formation. We have thus created series of initial conditions 
where one parameter only was varied, so as to best put for- 
ward the effect of this specific parameter. We considered in 
turn the effect of the scale- length of the mass distribution, of 
the cut-off radius of the galaxy distribution and of the mass 
fraction in the common halo. Most important, we have con- 
sidered both matter distributions with a core and matter 
distributions with a cusp, and have studied the effects of 
these two very different distributions on the evolution of the 
group and on the properties of the central galaxy. In other 
words we have been able to create our initial conditions such 
as to put forward best the phenomena we wished to study. 
Of course cosmological simulations and in particular obser- 
vations will be necessary to tell us whether all the spectrum 
of initial conditions considered here are possible or realistic. 
Our approach is similar to that used e.g. in studying the 
formation of bars in disc galaxies. It is possible to conclude 
that the growth rate of the bar is smaller in discs which 
are hot and/or have an important spheroidal component in 
their center or main parts, without having to worry about 
the exact form of the radial density profile of the disc or 
halo component. Once the essential physics of the bar for- 
mation is understood, then one can start worrying about 
whether this or that type of radial density profile is more 
realistic, follows best from cosmological initial conditions, 
or is in better agreement with the observations. 

The fraction of mass in the common halo is, in all our 
simulations, considerably smaller than what is believed to be 
the case in real groups or clusters. Nevertheless, the range 
explored is sufficient for us to establish how the fraction of 
mass in the common halo influences the evolution. It also 
allowed us to establish that the form of the radial profile of 
the central object in Plummer simulations is incompatible 
with observed profiles in cases with a high fraction of mass 
in a common halo. Simulations with a higher fraction of 
mass in the common halo would necessitate considerably 
more CPU time for two reasons. First the total number of 
particles in the simulations would increase considerably if 
the number of particles per galaxy is kept constant. Also 
the time it would take for the central object to form would 
also increase, thus leading to longer CPU times, if the time- 
step is kept constant. Nevertheless, it would be interesting 
to pursue in the future simulation with higher percentage 
of mass in the common halo, particularly in order to test 
the formation of the luminous halo in the outer parts of the 
central objects. 

Probably the most serious drawback of our simulations 
is the fact that our galaxy models are very simplistic. Thus 
we do not distinguish between luminous and dark matter, 
and have not considered any galaxies with a disc, with 
anisotropy or with rotation. Thus we can not distinguish 
between dark and luminous mass in the central object. Sim- 
ulations with more realistic galaxy models would necessi- 



tate a considerably higher number of particles, particularly 
if disc galaxies are considered. Such simulations would ne- 
vertheless be very interesting, and should be envisaged now 
that increases in computing power make them possible. 

We have also neglected the existence of gas, both cold 
and hot, both in the galaxies and in the common background 
halo. Thus effects like ram-pressure stripping of individual 
galaxies are totally neglected, but this should not infiuence 
the formation of the central object. It should of course be 
kept in mind that the fraction of common background halo 
in our simulations includes all mass encompassing the whole 
cluster and not linked to individual galaxies, i.e. it includes 
also the hot gaseous component. 

The above summarise the main simplifications and as- 
sumptions underlying our simulations. It is obvious that 
there is room for improvement, but none of our assumptions 
should alter the main physical results found in this paper. 



6 SUMMARY AND DISCUSSION 

We have performed N-body simulations of the dynamical 
evolution of groups of 50 galaxies embedded in a common 
background mass distribution. We considered two density 
distributions, the Plummer and the Hernquist one. The for- 
mer has a core, while the latter has a cusp. The main goal 
of this paper was to compare the evolution of the group and 
the properties of the central giant galaxies formed in these 
two very different matter distributions. 

As in the case of simulations of groups with no common 
background distribution of matter (GAG) , the evolution of 
the systems within a common halo is driven mainly by the 
merging instability (Carnevali et al. 1981). In all cases a gi- 
ant galaxy is formed in the central parts by merging of some 
secondary galaxies and by accreting material stripped from 
the rest of the galaxies by tidal forces. The evolution rate of 
the systems is very sensitive to the initial conditions and is 
controlled mainly by the density in the central parts. Dis- 
tributions with higher densities resulting either from higher 
central concentrations or higher densities over a larger area, 
have a faster evolution than distributions with lower densi- 
ties in the central parts. One should nevertheless remember 
that the dynamical time of denser configurations is shorter. 

To study more quantitatively the evolution of the dif- 
ferent simulations we have followed the time evolution of 
some global parameters, like the number of galaxies that 
have not yet merged, their mean distance from the center 
of the group, their mean mass and their velocity dispersion. 
The total number of galaxies decreases with time in all cases. 
This decrease is faster for configurations with a larger frac- 
tion of the total mass in galaxies. This temporal decrease is 
particularly important for the first 2 Gyrs of the evolution 
of cusped initial distributions, where there is a strong initial 
merging. In the Plummer simulations, where the initial dom- 
inant effect is the tidal stripping of material, the temporal 
decrease of the number of galaxies is only important after 
the first 2 Gyrs. The effect of the initial conditions on the 
mean distance of the galaxies from the center of the group, 
on their mean mass and on the velocity dispersion of the 
galaxy system is less spectacular. 

We have also studied the effect of substituting the live 
common background halo of particles by a rigid halo with 
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the same parameters by comparing two simulations with the 
same initial conditions, except that the one has a rigid and 
the other a live halo. The global evolution of the two simula- 
tions is quite different. In particular we find that the number 
of galaxies that have not yet merged decreases much faster 
in the simulation with the live halo. This can be easily un- 
derstood since the evolution in a live halo includes the effects 
of dynamical friction which brakes the galaxy motion and 
favours the merging with the giant galaxy. Thus, although 
simulations with rigid common halos require considerably 
less CPU, they should be avoided, since their results are not 
reliable. 

The different profiles of the mass distribution influence 
both the evolution of the group and the properties of the 
giant galaxy formed in the centre. In the case of a Plummer 
density distribution the individual galaxies undergo pairwise 
mergers giving rise to some larger galaxies in the core of 
the cluster. At some later step these galaxies merge and 
form the giant central galaxy. On the other hand, in the 
Hernquist case the high central concentration of galaxies in 
the center leads to a high '"burst" of mergers in the central 
parts during the initial steps of the evolution and the giant 
central galaxy is quickly formed. 

Not only the merging histories but also the quantity 
of stripped material from the satellite galaxies is different 
in the two cases. In the Plummer cases the fraction of the 
total mass in the central galaxy that comes from stripping is 
very important, especially early on in the simulation. In the 
Hernquist cases the fraction of the mass of the central galaxy 
that comes from stripped material is not so important, but 
amounts nevertheless to a considerable mass, since the total 
mass in the central object is important. 

These different evolutionary histories give rise to cen- 
tral objects with very different observable properties. We 
have studied the projected density distributions of the gi- 
ant galaxies at different time-steps and for different random 
projections. The shape of these profiles depends heavily on 
the initial distribution of the mass. In the cases of Plummer 
distributions with half the total mass in the common back- 
ground, the projected density profile of the giant galaxy does 
not follow the r^'** law at any time during the simulation. 
The radial density profiles of these central objects are not 
observed in the real giant galaxies in the center of the clus- 
ters. Only for the case where most of the mass is initially in 
the galaxies does the giant galaxy formed within a Plummer 
distribution follow the r^'*' law typical of elliptical galaxies. 

On the contrary, the galaxies formed within Hernquist 
distributions do follow r^'** laws in the main parts of their 
bodies in all the cases. Moreover, in the external parts, they 
show the typical projected density excess of cD galaxies in 
the case where a considerable fraction of the mass is initially 
in the galaxies. This excess over the r^'** law is not a pro- 
jection effect and thus reflects the real matter distribution 
of the central object. Furthermore, this excess gets larger as 
the dynamical evolution of the group proceeds. 

Thus our simulations suggest that central galaxies with 
r^'** profiles can form either if the fraction of mass in the 
common halo is very small, or if the density distribution 
is cusped. Since the amount of matter observed in the com- 
mon background is considerably larger than what would lead 
to r^''^ profiles for all density distributions, our simulations 



argue strongly in favour of cusped matter distributions, at 
least for those groups and clusters where a BCM has formed. 
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